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LUNAR TERRAIN MAPPING AND RELATIVE-ROUGHNESS 

ANALYSIS 


By Lawrence C. Rowan, John F. McCauley, and Esther A. Holm 


ABSTRACT 

Terrain maps of the equatorial zone (long 70° E.-70 0 W. 
and lat 10° N-10° S.) were prepared at scales of 1:2,000,000 
and 1:1,000,000 to classify lunar terrain with respect to 
roughness and to provide a basis for selecting sites for Sur- 
veyor and Apollo landings as well as for Ranger and Lunar 
Orbiter photographs. The techniques that were developed as a 
result of this effort can be applied to future planetary explora- 
tion. 

By using the best available earth-based observational data 
and photographs 1: 1,000,000-scale and U.S. Geological Sur- 
vey lunar geologic maps and U.S. Air Force Aeronautical 
Chart and Information Center LAC charts, lunar terrain was 
described by qualitative and quantitative methods and divided 
into four fundamental classes: maria, terrae, craters, and 
linear features. Some 35 subdivisions were defined and map- 
ped throughout the equatorial zone, and, in addition, most of 
the map units were illustrated by photographs. 

The terrain types were analyzed quantitatively to charac- 
terize and order their relative-roughness characteristics. 
Approximately 150,000 east-west slope measurements made 
by a photometric technique (photoclinometry) in 51 sample 
areas indicate that algebraic slope-frequency distributions are 
Gaussian, and so arithmetic means and standard deviations 
accurately describe the distribution functions. The algebraic 
slope-component frequency distributions are particularly use- 
ful for rapidly determining relative roughness of terrain. 

The statistical parameters that best describe relative 
roughness are the absolute arithmetic mean, the algebraic 
standard deviation, and the percentage of slope reversal. 
Statistically derived relative-relief parameters are desirable 
supplementary measures of relative roughness in the terrae. 
Extrapolation of relative roughness for the maria was demon- 
strated using Ranger VII slope-component data and regional 
maria slope data, as well as the data reported here. It appears 
that, for some morphologically homogeneous mare areas, 
relative roughness can be extrapolated to the large scales from 
measurements at small scales. 

INTRODUCTION 

SCOPE AND PURPOSE 

In October 1963, the U.S. Geological Survey be- 
gan, on behalf of the National Aeronautics and 
Space Administration (contract No. W-11,775), a 


program of terrain studies of the lunar equatorial 
zone (lat 10° N.-10° S., long 70° E.-70 0 W.), the 
area most suitable from an engineering standpoint 
for early Apollo manned landings (fig. 1). These 
studies complemented the 1 :1, 000, 000-scale lunar 
geological mapping program and aided in selection 
of targets for Ranger VII and VIII, Lunar Orbiter, 
and Surveyor I missions. 

The principal objectives of the terrain analysis 
program were (1) to examine the roughness charac- 
teristics of the existing lunar geologic map units, 
(2) to redefine and regroup these geologic units into 
mappable terrain or morphologic units without re- 
gard to their stratigraphic position, (3) to charac- 
terize the relative roughness of these terrain units 
by means of simple statistical expressions and to de- 
pict their distribution in a series of terrain maps 
which could be used for mission planning, (4) to 
study the smoothest parts of the lunar equatorial 
belt visually at the telescope to evaluate their rela- 
tive suitability as landing areas for Surveyor and 
Apollo spacecraft and thereby further narrow the 
area to be photographed by Lunar Orbiter space- 
craft, and (5) to refine and automate existing photo- 
metric slope-measuring techniques for the purpose 
of generating quantitative slope data for statistical 
analysis. 

Many reports were prepared to meet flight sched- 
ules of the unmanned Lunar Orbiter and Surveyor 
programs; most have not been formally published. 
Some of these, released in open files by the U.S. Geo- 
logical Survey, are listed under “Selected refer- 
ences.” The integration of these preflight support 
studies into the mission plans of the early Surveyor 
and Orbiter experiments was an element in the 
overall success of the missions. 

The purpose of this paper is (1) to summarize 
lunar terrain mapping techniques described in these 
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Figure 1. — U.S. Air Force Chart and Information Center photomosaic showing the extent of the lunar terrain study area. 


mostly unpublished sources, (2) to describe the pho- 
tometric method used to extract slope data from tel- 
escopic photographs, (3) to describe the various 
statistical parameters that were applied to the slope 
data to establish relative roughness, and (4) to re- 
view, briefly, the major applications to targeting of 
the unmanned lunar probes, particularly Surveyor I 
and the early Lunar Orbiters. Since completion of 
this work, the high-resolution data from the many 
successful probes have rendered most of the tele- 


scopic measurements and maps obsolete. For this 
reason emphasis is placed on the techniques and 
methodology developed through mid-1966, time of 
the first Lunar Orb iter flight, rather than on specific 
results. The techniques remain applicable to higher 
resolution terrain analyses of the moon and other 
terrestrial planets. The slope measurements and sta- 
tistical classifications are currently in various stages 
of refinement and engineering application because 
of a continuously improving stream of data. 
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PREVIOUS WORK AND SOURCES OF DATA 

An antecedent of the program described here was 
the pioneering work of Hackman and Mason (1961) . 
Their report consists of the following: (1) A de- 
scription of physiographic divisions of the moon, 
which divides surfaces features according to “ex- 
tent of preservation, type of modification, type of 
surface material, elevations, slopes and structural 
disturbances/' (2) a generalized photogeologic map 
of the moon, which classifies surface materials into 
premaria, maria, and postmaria ages, (3) a map 
showing the distribution of lunar rays, which are 
interpreted as ballistic ejecta from the larger impact 
craters, and (4) a brief interpretative text and a 
unit-by-unit discussion of the predicted landing, 
mobility, and construction characteristics within 
each of the physiographic divisions. The terrain 
study described here is thus an extension of Hack- 
man and Mason's work and represents a similar ap- 
plication of geologic methodology to the engineering 
aspects of planetary exploration, particularly the 
problem of landing-site selection. Additional terrain 
studies directed more towards surface mobility prob- 
lems have since been undertaken by Rozema (1968) . 

Much of the topographic and geologic data utilized 
in the present study was produced by two systematic 
1:1, 000, 000-scale telescopic mapping programs: (1) 
the U.S. Geological Survey lunar geologic mapping 
program, summarized by Wilhelms (1970) and (2) 
the 1:1, 000, 000-scale topographic mapping program 
of the U.S. Air Force Aeronautical Chart and Infor- 
mation Center. Additional sources of raw data in- 
cluded (1) published and unpublished lunar photo- 
graphs from Mount Wilson, Lick, McDonald, Pic du 
Midi, and Catalina Observatories and (2) visual 
studies of selected parts of the lunar equatorial belt 
at the U.S. Geological Survey 30-inch reflecting tele- 
scope near Flagstaff, Ariz.; the Kitt Peak National 
Observatory 60-inch reflecting telescope and the Mc- 
Math solar telescope, both west of Tucson, Ariz.; 
and the Lick Observatory 36-inch refracting tele- 
scope, Mount Hamilton, Calif. 

Quantitative topographic data consisted of (1) 
east-west slope-component measurements obtained 
by photometric techniques developed by the U.S. Geo- 
logical Survey (Wilhelms, 1964; McCauley, 1964; 
Rowan and McCauley, 1966) and (2) relative-relief 
measurements made by the shadow technique (Ko- 
pal, 1962) and furnished by the U.S. Air Force Aero- 
nautical Chart and Information Center, St. Louis, 

Mo. 
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LUNAR TERRAIN MAPPING 

Terrain maps generally portray the areal distribu- 
tion patterns of relatively homogeneous land-surface 
forms, such as plains, tablelands, hills, and moun- 
tains, without genetic implication. Examples of this 
type of map, adapted from the work of Hammond 
(1964), are sheets 61-64 of the National Atlas (U.S. 
Geol. Survey, 1968). Large-scale terrain maps, such 
as those prepared by the Department of the Army 
(Chief of Engineers, 1958), usually serve a specific 
practical purpose, such as trafficability analysis. The 
terrain, or land-surface form, map is distinct from 
the physiographic map (for example, sheet 59, Natl. 
Atlas), which attempts to provide an understanding 
of landform genesis, and from geologic and tectonic 
maps, wherein the fundamental classification is ac- 
cording to geologic age or the role of a body of rocks 
in the tectonic history of a region. 

The lunar terrain maps prepared during this study 
are special-purpose landform maps that depict vary- 
ing degrees of suitability for both unmanned and 
manned landings. They convey no stratigraphic or 
genetic information ; however, because the Geo- 
logical Survey's 1:1, 000, 000-scale lunar geologic map 
series was used extensively in the compilation of the 
terrain maps, the terrain unit boundaries generally 
coincide with stratigraphic unit boundaries, particu- 
larly where a given stratigraphic unit is character- 
ized by unique terrain characteristics. Generally, 
terrain unit boundaries depart from stratigraphic 
contacts where morphologic expression on two or 
more stratigraphic units is similar or where the 
morphologic properties of a stratigraphic unit vary 
laterally, as around the younger multiring basins. 
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EXPLANATION 


Median slope frequency ranges from zero (no resolvable relief) for small areas to 1.5° as determined by 
photometric techniques. Generally consists of the darker parts of the maria 


Median slope frequency ranges from zero (no resolvable relief) for small areas to 1.5 as determined by photo- 
metric techniques. Generally consists of the lighter parts of the maria and ray-covered mare areas. Visual 
studies indicate that these areas are, under the conditions of best resolution, rougher than the darker parts of 
the maria 


Median slope frequency ranges from 1° to 2°, as determined by photometric techniques. Consists of mare 
areas that show numerous structural elements such as sinuous ridges, domes, rilles.or chain craters which 
contribute to an undulating, uneven topography 


Median slope frequency, as determined by photometric techniques, ranges from 1.5° to 2.5°. Consists generally 
of relatively smooth basins occupying the floors of the largest craters 


Median slope frequency, as determined from shadow progression and photometric techniques, ranges from 2° 
to 6°. Consists of individual craters superimposed on smoother, surfaces. Includes both the crater floor, 
walls, and rim deposits to a distance of approximately one crater radius 


Median slope frequency, as determined from shadow progression, and photometric techniques, ranges from 3° 
to a measured maximum of 8°. Surface consists primarily of a complex sequence of overlapping craters and 
crater deposits; the mean slope value of the interior walls of individual crater is 37 


Figure 2.— Terrain of part of the Sinus Medii area, mapped at 1:2,000,000 scale; units described by slope-component statistics (after McCauley, 

1964, p. 27-30). 


CONTRIBUTIONS TO ASTROGEOLOGY 








LUNAR TERRAIN MAPPING AND RELATIVE-ROUGHNESS ANALYSIS 


G5 


1 :2, 000,000-SCALE TERRAIN MAPPING 

The first effort at lunar terrain classification re- 
sulted in a six-unit 1:2, 000, 000-scale map of the area 
from long 15° E.-60 0 W., and from lat 10° N.-10 0 S. 
(McCauley, 1964). The central part of this map is 
reproduced in figure 2. The six terrain units were 
defined and ranked on the basis of median slope val- 
ues (McCauley, 1965), which were calculated from 
some 40,000 measurements made for slope elements 
about 0.75 km in length. In order of decreasing 
roughness, the units are as follows : ( 1 ) Rugged terra, 
consisting mostly of closely spaced hills and ridges, 
degraded craters, and crater remnants with collective 
median slope values from 3° to as much as 8°. (2) 
Rim, wall, and floor deposits of large fresh craters. 
Median slope values in this unit range from 2° to 6°, 
but slopes on interior walls may be as much as 30°. 
(3) The relatively smooth-appearing floors of the 
large circular basins in the terrae. These materials 
are higher in albedo than most of the material in the 
more extensive mare basins, but are similar other- 
wise in their telescopic properties; median slope val- 
ues range from about 1.5° to 2.5°. (4) Ridges, domes, 
and rilles in the maria; these areas have median slope 
values from about 1° to 2°. (5) Ray-covered mare 
surfaces; these areas typically surround large fresh- 
appearing craters such as Kepler and Copernicus and 
can be traced radially from the parent crater for 
distances of 10-20 crater diameters. They do not 
exhibit, at average telescopic resolution, any greater 
telescopically measurable roughness than the non- 
rayed parts of the maria. Visual studies, however, 
indicate that at very low angles of illumination the 
ray-covered mare areas appear to be somewhat 
rougher than nonrayed mare areas. (6) Nonrayed 
maria; these areas generally show little or no resolv- 
able relief, and median slope values are consistently 
less than 1.5° in the intercrater areas. 

Prior to this terrain categorization it was not clear 
to most lunar exploration planners that the rayed 
parts of the maria are rougher than the nonrayed 
areas or that the floors of the terra basins are meas- 
urably rougher than the maria. Demonstration of the 
greater roughness of these units led to the general 
realization that the nonrayed parts of the maria 
should be examined first by unmanned spacecraft. 

1:1, 000,000-SCALE TERRAIN MAPPING 

During compilation of the 1:2, 000, 000-scale terrain 
map, it became apparent that a more detailed map 
of the equatorial zone would be needed for both Sur- 
veyor and Lunar Orbiter site-selection studies, par- 
ticularly in complex areas. Compilation of a more 
detailed map at the 1:1,000,000 scale became feasible 


with the later acquisition of some 65 photographs 
taken in 1964 at the 120-inch telescope at Lick Ob- 
servatory. The resolution of these photographs is 
generally better by a factor of two than previously 
available photographs (Kuiper, 1960, 1961; Whitaker 
and others, 1963). Continued refinements in available 
geologic maps and topographic charts also facilitated 
1:1, 000, 000-scale terrain mapping. 

Figure 3 shows part of one of the six 1:1, 000, 000- 
scale sheets that cover the area from lat 10° N. to 10° 
S. and from long 60° E. to 60° W. Table 1 gives de- 
scriptions for each unit, and figures 22-29 show tele- 
scopic photographs of most of the terrain units. The 
1:1, 000, 000-scale terrain Glassification is based on the 
recognition of four fundamentally different types of 
surface forms: (1) mare, (2) terra, (3) craters, and 
(4) linear features such as ridges, rilles, rays, es- 
carpments, and mare domes. These four morphologic 
types are in turn divided into related variants. 

Table 1 . — Terrain units of the lunar equatorial belt 

I— MARE: 

I-A Mare, smooth plains mostly in and around depressed 
parts of large circular basins. 

I- B Mare, dark, smooth plains mostly near margins of large 
circular basins and along west edge of Oceanus Procel- 
larum. 

II— TERRA: 

II- A Terra, smooth plains of intermediate albedo with 
abundant ^l>km craters in terra basins. 

II-B Terra, gently rolling to hummocky areas. 

II-C Terra, hummocky to moderately blocky areas. 

II- D Terra, rough, high blocky areas with deep intervening 
linear troughs; mostly marginal to south rim of Mare 
Imbrium. 

III— CRATERS: 

III- A Topographically sharp well-formed craters with or 
without associated bright rays : 

III-A-1 Wall (interior), undifferentiated: 

III-A-la Steeply sloping, rough terraced or scalloped 
walls or detached blocks of wall material. 
III-A-lb* Wall, moderately smooth but steeply slop- 
ing. 

III-A-2 Floor, undifferentiated : 

III-A-2a* Floor, smooth plains with little local relief, 
generally of limited extent. ■ 

III-A-2b Floor, hummocky to hilly, rolling topog- 
raphy. 

III-A-2c Peaks, generally at or near the center of 
the floor; rough steeply sloping flanks. 

III-A-3 Rim, undifferentiated: 

III-A-3a* Rim, rough, hummocky or with concentric 
troughs and ridges, generally proximal to the rim 
crest. 

III-A-3b* Rim, topography consisting .of radial 
ridges and troughs, distal to hummocky area on 
rim. 

III-B Modified craters, generally of subdued form: 


♦See footnote at end of table. 
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Table 1. — Terrain units of the lunar equatorial belt — Cont. 

III — CRATERS — Continued 

III-B-1 Wall (interior), undifferentiated, but generally 
smooth. 

III-B-2 Floor, undifferentiated: 

III-B-2a* Floor, smooth areas generally within large 
markedly deformed craters on terra. 

-• III-B-2b* Floor, rough; marked by blocks and irreg- 
ular terraces; hummocky. 

III-B-2c* Peaks, centrally located; height above crater 
floor is generally less than that of III-A-2c peaks. 

III- B-3* Rim, subdued rolling to hummocky topography; 
around small craters, this topography cannot be fully 
resolved. 

III- C Crater fields and craters 2-10 km in diameter. 

IV — LINEAR FEATURES : 

IV- A Ridges, both linear and sinuous, located within the 
mare. 

IV-B Domes, some have small summit craters barely resolv- 
able at 1-km resolution; some domes appear rough tex- 
tured at the limit of telescopic resolution. 

IV-G Rilles (rimae) : 

IV- C-1 Rilles, linear, occur in both mare and terra ; gen- 
erally^ have sharply defined steeply sloping interior 
walls. 

v- * TV-C-2* Rilles, sinuous; generally occur in terra but also 
along mare margins ; commonly originate in crater ; 
some terminate in Small bifurcating rilles. 

^IV-D Plateaus,: smooth to gently, rolling surfaces, bounded 
by escarpments; generally within mare. 

IV-E* Escarpments (rupes). 

IV— F*. Troughs. 

IV-G Depressions, irregular or with ill-defined boundaries. 

IV-H Chain craters. 

♦Units, not shown in figure 3. 

In general, the maria are lower in elevation than 
the terrae and also have surfaces of lower relief and 
reflectivity (albedo) . Wide troughs and plateaus, bor- 
dered by scarps in some places, together with ridges, 
rilles, domes,, and craters lend relief to and disrupt 
vthe' topographic continuity of these dark plains. If 
the linear features and craters were removed, the 
surface would be almost completely featureless at 
typical telescopic resolution (1 km). The most wide- 
spread mare type is I-A (table 1). The albedo of the 
maria is locally, however, very low, and the term 
‘‘mare, dark” (table i, I-B) has been used to desig- 
nate such areas (fig. 3), which generally have dem- 
onstrably fewer small craters. Relatively smooth 
plains also occur in the terrae within circular basins 
or degraded large craters; these plains differ from 
mare plains by having a higher density of craters 
that are at and just below the limit of telescopic 
photography, but visible on Lunar Orbiter photo- 
graphs. In the Sinus Medii area, depicted in figure 3, 
most of the mare is of type I-A, with a density of 


craters and mare ridges that is reasonably repre- 
sentative of other maria areas on the earthside 
hemisphere. 

The terrae in the central part of the equatorial belt 
are morphologically complex, even at telescopic reso- 
lution. Nevertheless, a fourfold division (II-A, II-B, 
II-G, II-D) evolved after several earlier, unsuccess- 
ful classification attempts. Lunar Orbiter photo- 
graphs show that the terrae are even more complex 
than had been realized, but the original four divisions 
are sufficiently general to remain applicable within 
the map area. 

Terra plains of intermediate albedo, II-A, are 
common at the mare-terra boundaries and also in 
basins of structural origin or in floors of highly de- 
graded craters (fig. 3). Hummocky to faintly braided 
terrain, II-B, is widely distributed around the south 
boundaries of Mare Imbrium and is quite different in 
texture from the other units. Unit II-C is character- 
ized by subdued linear ridges and troughs (fig. 3) 
oriented radially to Mare Imbrium. The roughest 
terrain unit, II-D, consists of strongly lineated terra 
with isolated positive blocks, steep ridges, and 
troughs. 

Craters, terrain class III, are divided into two 
main categories: III-A, topographically sharp well- 
formed craters with or without associated bright 
rays and III-B, subdued more irregular crater forms, 
including most of the larger circular features in the 
terrae. In each of these categories, morphologic sub- 
units can generally be distinguished within larger 
craters ( > 10 km). However, the smaller craters, 
2-10 km in diameter, cannot be classified as to type 
and are, therefore, placed in an undivided unit, III-C 
(fig. 3). Large relatively unmodified craters, such as 
Copernicus and Theophilus, exhibit a wide 1 range of 
roughness characteristics at 1 km resolution, depend- 
ing on location within a subunit. The interior walls 
are clearly the least hospitable terrain unit, with 
slopes of as much as 30°. 

Linear features, unit IV, constitute the fourth 
group of major terrain units. Ridges, rilles (both 
straight and sinuous), escarpments, chain craters, 
troughs and elongate plateau areas, and faults and 


Figure 3. — Terrain of a part of the Sinus Medii area mapped 
at 1:1,000,000 scale (after Holm and others, 1965). Compare 
with figure 2. Table 1 gives descriptions of each unit. Light 
dashed lines indicate approximate unit boundaries; heavy 
dashed lines locate narrow linear depressions. Heavy, solid 
lines indicate faults; solid lines with diamonds and arrows 
locate ridge crests. Stippled pattern locates areas of ray 
deposits. 
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lineaments each have individual and identifiable top- 
ographic expression. The last two are not mappable, 
however, as units at this scale because of their lim- 
ited areal extent. (See types of terrain units, figs. 
22-29.) The frequency of any one of these types 
within an area may be such that it dominates the 
terrain pattern and thus markedly alters the surface 
roughness. Also included, mostly for convenience, in 
this category are the round to elliptical domes pres- 
ent locally in the maria. 

VISUAL TELESCOPIC TERRAIN MAPPING 

Craters as small as a half a kilometer and linear 
clefts as narrow as 100 meters can be seen with a 
medium-aperture telescope when viewing at low illu- 


mination angles. This resolution represents nearly a 
twofold improvement over the terrain detail on most 
high-quality lunar telescopic photographs. For this 
reason a program of visual terrain studies of small 
selected areas was conducted using primarily the 
Lick Observatory 36-inch telescope. Ten mare sites 
were picked on the basis of their apparent smooth- 
ness on existing telescopic photographs. Figure 4 
shows a map of one of these study areas, north of 
the crater Flamsteed. It is representative of the vis- 
ual mapping approach and is of additional interest 
because it was utilized in the final selection of the 
Surveyor I landing site. The terrain of the Flamsteed 
area was divided by. D. F. Crowder (in McCauley, 
1964, p. 36) into five units arranged in order of rela- 



EXPLANATION 


Craters, including the rim deposits for a dis- 
tance of one crater radius. Interior slopes 
average 37° 



Broad ridges with slopes from 1° to 3° 
Roughness appears greatest along crests 
and in lower parts of structures 


E3 

Relatively flat but surficially rough material, 
as determined from the percentage of 
shadowed area at low sun illuminations 



Moderately smooth mare material with few 
alternating lighted and shadowed areas 



Smoothest terrain corresponding to the 
darkest part of the maria in the region 

□ 

Unstudied areas 


Contact between terrain units 
Dashed where indefinite 


Figure 4. Sample lunar terrain map prepared by visual observation at the Lick 36-inch refracting tele- 
scope. (After McCauley, 1964, fig. 10.) 
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tive roughness. Roughness was estimated by the 
incidence of small resolvable shadows close to the 
terminator. Regions showing the largest percentage 
of shadowed terrain alternating with illuminated 
areas were judged to be the roughest, and those 
showing the least variation in illumination were con- 
sidered the smoothest. The smoothest areas are gen- 
erally coincident with dark mare — a relationship 
found to be true in the other areas studied visually. 
Focus on these darker mare areas during the site- 
selection phase of the Lunar Orbiter program facili- 
tated accomplishment of its primary mission by the 
third flight and freed the two remaining spacecraft 
for moon-wide reconnaissance and detailed examina- 
tion of areas of scientific interest. 

One of the major results of the visual terrain 
studies program was the discovery that ray-covered 
mare areas are rougher than nonrayed mare areas. 

SLOPE-MEASURING PROGRAM 

Strahler (1950, 1954, 1956) demonstrated that 
slope is the most important element in quantitative 
landform classification; Wood and Snell (1959) also 
showed that slope data are more effective than other 
topographic quantities for the purpose of terrain 
classification. Numerous other measures of terrain 
geometry have been used in terrestrial studies to 
completely describe the topography, but slope sta- 
tistics seem most appropriate for relative-roughness 
studies. 

Methods for obtaining reasonably accurate lunar 
topographic data have existed since the telescopic 
work of William Herschel and, later, that of Schroter 
in the 18th century (Kopal, 1962). The total quantity 
of reliable topographic data at the outset of the ter- 
rain studies program was, however, quite meager. 
The only significant body of useful topographic in- 
formation consisted of the relative-relief measure- 
ments made by the U.S. Air Force Chart and Infor- 
mation Center, St. Louis, Mo., in support of their 
1:1, 000, 000-scale lunar topographic mapping pro- 
gram. 

One of the first tasks of this program then was to 
improve existing slope-measuring techniques so that 
the major recognizable terrain units could be studied 
quantitatively. The photometric technique first de- 
scribed by van Diggelen (1951) was selected as the 
most promising from the standpoint of providing a 
large quantity of slope data in the shortest possible 
time. It also has the advantage of being most sensi- 
tive in areas of low relief and uniform albedo, the 
type of terrain most suitable for spacecraft landings. 
After considerable refinement, the van Diggelen tech- 


nique was used systematically to generate about 
150,000 individual slope measurements, which were 
later statistically analyzed and used for ordering 
major terrain units according to their relative 
roughness (Rowan and McCauley, 1966). 

The term “photoclinometry” has been used widely 
in recent years to describe this method of deriving 
topographic slope information from lunar imagery 
(photographic or electronic). The term distinguishes 
between the techniques of photogrammetry, on the 
one hand, and the application of lunar photometry 
to the derivation of topographic data, on the other. 
As used herein, it is restricted to that aspect of lunar 
photometry concerned directly with the extraction 
of slope data from monoscopic pictures. Watson 
(1968) described the technique in detail, emphasizing 
its use in extracting slope data from spacecraft 
images. 


PHOTOCLINOMETRY 

The photoclinometric technique depends on a pre- 
cise knowledge of the lunar photometric function, 
which is an expression of the reflectance properties 
of the lunar surface. The geometric conditions for 
earth-based photoclinometry are considerably less 
complex than those obtained when viewing the lunar 
surface from spacecraft, as described by Watson 
(1968). The description given here will be limited to 
the simpler case of earth-based observation. 

The lunar surface is a strong back-scatterer of in- 
cident visible radiation (Minnaert, 1961, p. 222-233), 
and reflectance is greatest in the direction of the 
incident beam. The brightness at any point on the 
lunar surface is a function of the normal albedo and 
the geometry of the incident and emergent light, so 
B = f(i, e, g)n 

where 

B = brightness; 

i = the angle of incidence, the angle between 
the incident ray from the sun and a line 
perpendicular to the local surface; 
e = the angle of emergence, the angle between 
the emergent ray to the sensor and a line 
perpendicular to the local surface; 
g = the phase angle, the angle between the inci- 
dent ray and the emergent ray (that is, 
the sum of Li and Le), measured in a 
plane containing both rays — the phase 
plane; 

/ ( i , e, g) = the photometric function; and 
n = normal albedo. 

These relationships are shown in figure 5. 

Lunar photographs taken at oblique angles of illu- 
mination show a close correspondence between 
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Figure 5. — Geometry of the lunar photometric function, 
where i — angle of incidence ; e = angle of emergence ; and 
g = phase angle in the special case where all lie in the 
same plane. 

selenographic longitude (or sun elevation) and iso- 
photes (lines of equal brightness) near the equator 
(Minnaert, 1961, p. 222-233). This correspondence 
permits the measurement of slope components in an 
east-west direction, but provides no slope informa- 
tion in directions not included in the phase plane. 
Figure 6, a polar view of the moon, illustrates how 



Figure 6. — Polar view of the moon, showing the geo- 
metric relations between selenographic longitude and 
east-west slope components. (After van Diggelen, 
1951.) 

the slope of a surface element in an area of uniform 
albedo can be calculated ; f or areas where the illumi- 
nation conditions are&nowh.-. Consider two topo- 



Figure 7. — Idealized microphotometer traverse on the lunar 
surface and the standard curves for different relative 
albedo units for a particular lunar photographic plate 
(negative). The excursions on the traverse are due to 
both albedo and local slope effects. (After McCauley, 1964, 
fig. 3.) 

graphic elements, P and Q, of similar albedo and 
equal brightness because of their similar angular 
relations to the sun; element Q is horizontal or tan- 
gent to the curvature of the surface, where element 
P is sloping. The angle of slope of element P, LAEB 
(measured from a tangent to the surface), is given 
by LCOQ, the difference in selenographic longitude 
between the two elements. 

LAEB = LCPD = LCOQ 

Brightness changes in the east-west direction re- 
sulting from changes in the angle of slope can, there- 
fore, be directly equated with changes in sun angle 
or selenographic longitude, LCOQ, for topographic 
elements near the equator. Slope components, in the 
east-west direction can, therefore, be measured if the 
photometric function and the effects of albedo varia- 
tion from element to element can be eliminated. 
Brightness values that deviate from the average pho- 
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tometric function are the result of local slope either 
toward or away from the sun. The elements sloping 
toward the sun are brighter for a particular sun 
elevation angle than elements that are horizontal; 
conversely, elements sloping away from the sun are 
darker. The increase in brightness is a measure of 
the degree of slope determined by comparing the 
measured value with an empirically derived photo- 
metric curve for the particular photograph. 

Negative glass copy plates were used on a digitized 
Jarrel-Ash microphotometer to acquire the lunar 
slope data. Since the intensity on a negative lunar 
plate increases exponentially toward the terminator 
(fig. 7), a semilogarithmic plot of intensity versus 
sun elevation (or distance from the terminator in 
selenographic longitude) can be prepared for each 
plate by systematic sampling of areas determined by 
inspection to be smooth and of the same albedo. The 
resulting curve combines the effects of the lunar pho- 
tometric function and the Herter and Dryfield curve 
(H and D) of the plate itself. A number of curves, 
one for each recognizable albedo unit, were derived 
from the measured intensity values by a least-squares 
fit; they are referred to herein as standard curves. 
Departures in intensity from these curves indicate 
local slope effects in an east-west direction. Figure 7 
shows how an individual reading, X, within relative 
albedo unit 3, is compared with a point of equal in- 
tensity, Y, on the standard curve, to determine the 
slope of X, which is 2.2°. Large-intensity changes 
such as those between points A and B, and C and D 
(fig. 7) are typically the result of albedo differences, 
whereas changes such as those between points E and 
F, and F and G generally reflect changes in slope. 
Individual intensity readings are referred to their 
appropriate standard curve for slope determination 
by means of coding introduced during the plate 
traverse by the microphotometer operator. The inten- 
sity readings, referred to the standard curve for the 
particular albedo unit, are then processed by com- 
puter to acquire the actual slope-component data. 

Slope determination in areas of varying albedo re- 
quires prior categorization of the sample area into 
relative albedo units. A photographic technique de- 
veloped by Hawkins and Munsey (1963) and applied 
by McCauley (1964) to the problem of objectively 
discriminating relative albedo units was used. The 
technique consists of successively exposing a full- 
moon photograph to a very high contrast copy film 
(Kodak Ortho-Litho S). The black and clear response 
of the film to different exposure is utilized to outline 
isopleths of relative albedo, which then permit the 
derivation of the previously described standard curve 


within each albedo unit outlined. 

A detailed description of the methodology, the in- 
struments used, the logic of the computer programs 
employed, and some of the general results of the U.S. 
Geological Survey slope-measuring program were as- 
sembled by D. A. Bremmer (written commun., 1964). 
This work, along with a description of the first exten- 
sion of the method to nonmare areas by Wilhelms 
(1964) is a useful supplement to the published 
Reference already cited. 

LIMITATIONS OF EARTH-BASED PHOTOCLINOMETRY 

Although no systematic error analysis was made, 
limited comparisons of individual slopes measured 
photoelinometrically and by use of shadow-length 
measurements agree to within 10 percent in terrain 
sloping between 3° and 10°. In addition, profiles com- 
piled from photoclinometric measurements agree 
with those prepared from available U.S. Air Force 
Chart and Information Center LAC charts. 

The principal sources of error in earth-based 
photoclinimetry are (1) non uniformity of expo- 
sure and of development of the photographic plate, 
(2) variation in resolution within the plate 
caused by astronomical “seeing effects,” (3) errors 
in determination of local albedo variation which may 
be falsely interpreted as slope effects, (4) instrumen- 
tal errors in the microphotometer readings, 1 (5) er- 
rors in selection of the level areas that determine the 
standard curve for each albedo unit,, and (6) selec- 
tive censoring of the slope data in shadowed areas. 
The only measurement possible in shadowed areas is 
that value indicating the slope angle is greater than 
the sun angle at that point. Thus, the higher' slope 
values along any profile tend to drop out of the slope 
population as the terminator is approached. 

Instrumental errors . can be controlled ;, however, 
the relative contributions of each of the other sources 
of error to individual slope-component measurements 
cannot be accurately assessed. Clearly the most se- 
rious of possible errors stems from the lack of suit- 
able albedo correction. Figure 7 illustrates this prob- 
lem. If the reading at “X” were referred improperly 
to albedo unit 2 rather than 3, a 100 percent error in 
slope measurement would result. Additional unavoid- 
able sampling problems that contribute ambiguity to 
individual measurements are discussed later in this 
report. 

These limitations are not as serious as they first 


1 One such problem results from the slightly lesser measuring capability of 
the microphotometer in areas of high density than in areas of low density on 
the photographic plate; thus, accuracy of slope-component determination is 
slightly lower for areas sloping toward the sun than away from it on the 
negative plates used. 
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appear if the main purpose of the measurements is 
kept in mind. The measurements are not used to pro- 
duce topographic maps, but serve only as the raw 
material necessary for statistical evaluation of rela- 
tive roughness. Erroneous individual slope values 
are not as important as they would be if topographic 
mapping were the objective. All the aforementioned 
types of error can affect slope measurements in each 
part of the moon sampled, and the overall cumulative 
error for any block of measurements tends to be uni- 
form. Careful plate selection, determination of albedo 
variation, microphotometer calibration, and deter- 
mination of standard curves can keep the total error 
generally constant so that the relative order of meas- 
urable roughness of the units sampled is not changed. 
Statistically speaking, then, the lunar surface can 
be described and arranged into progressively rougher 
units, the relative order of which will not be affected 
by systematic slope-measuring errors. Owing to these 
errors, the final measurements may be consistently 
too high or tod low, but the order in which units are 
ranked will be unaffected. Preliminary comparisons 
of statistical slope values computed from profiles de- 
rived (1) photogrammetrically by S. S. C. Wu and 
(2) photoclinometrically by Baerbel Lucchitta indi- 
cate that the photoclinometric slope measurements 
are 15-25 percent lower than the photogrammetric 
ones. The profiles compared were derived, from 
Apollo 8 photographs and will appear in the pub- 
lished mission report. 

RELATIVE RELIEF 

The photoclinometric technique is relatively inef- 
fective for measuring slopes of rough terrain in the 
terminator region because horizontal surfaces are 
rare and large areas are in shadow; these factors 
preclude accurate derivation of the standard -curves 
and causes censoring of the steeper slope values. For 
these reasons, relative-relief measurements were 
chosen for statistical analysis of the rougher parts 
of the equatorial zone. 

Relative relief is the vertical displacement of a 
topographic feature, as measured along a horizontal 
line. The 1,200 measurements used here were ob- 
tained from William Cannell of the U.S. Air Force 
Chart and Information Center, St. Louis, Mo.; they 
were made by the shadow technique developed by 
Kopal (1962). Unlike the slope-component measure- 
ments, relative-relief measurements are not confined 
to 0.75 km slope lengths, but range from about 3 to 
5 km. They are, therefore, not strictly comparable, 
in terms of resolution, with the values obtained 
photoclinometrically. 
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Figure 8. — Locations of the 61 slope-component sample areas in the lunar equatorial zone. (After Rowan and McCauley, 1966, fig. 6.) 
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STATISTICAL ANALYSIS 

SAMPLING PROCEDURE 

The engineering applications of ranking terrain 
according to relative roughness for landing-site selec- 
tion studies dictate that only standard and easily 
computed statistical measures be examined. The use- 
fulness of these measures for terrain ranking de- 
pends on the accuracy of the measurements them- 
selves and on the degree to which the samples repre- 
sent the real slope populations at the scale of interest. 
The estimated accuracy of the measurements has 
already been discussed, and this section will deal 
mainly with the sampling procedures, the inherent 
problems therein, and the effects of these problems 
on relative-roughness rankings. 

Consider two general scale-dependent sample re- 
gimes: the terrain map scale, on which the smallest 
classified features are on the order of tens of kilom- 
eters, and the scale of the individual slope-component 
measurements, each of which is approximately 0.75 
km long. At the 1:1,000,000 map scale, terrain units 
are delineated visually on the basis of common suites 
of morphologic characteristics; thus, these units are, 
at telescopic resolutions, nearly morphologically ho- 
mogeneous and, therefore, have characteristic slope 
distributions and relative-relief populations. In most 
samples, the percentage of area occupied by units 
other than the populations of primary interest is 
small, so the slope-component statistics are reason- 
ably representative of the primary population. Some 
of the 51 slope sample areas (fig. 8) were, however, 
distinctly heterogeneous and thus had to be divided 
to obtain morphologic homogeneity. 

The most significant sampling problems arise from 
limitations of the photoclinometric technique and of 
the measuring equipment used. The photoclinometric 
technique, as previously shown, is restricted to east- 
west measurements, so true slope (maximum slope) 
measurements are made only when the trend of strike 
of the sloping surface is north-south; measurements 
of slopes that do not trend north-south represent a 
component of the true slope. Unidirectional sampling 
of this type would have little or no effect if the lunar 
surface were geometrically random, but many topo- 
graphic features of the lunar surface are preferen- 
tially oriented. Where the trends of these features 
are known, corrections for the departure from true 
slope can be made (van Diggelen, 1951), but were 
not because of the increased computational burden 
and the judgment that such corrections are not nec- 
essary for the type of statistical ranking attempted. 
Thus, most of the measurements used here are slope 
components, rather than true slope measurements, 


and the calculated statistics used for ranking rough- 
ness are systematically biased toward lower values. 
The bias, however, appears to be generally uniform 
and, therefore, does not seriously detract from the 
validity of the statistical terrain rankings. 

The sample length and the geometry of the terrain 
also determine the degree to which samples represent 
real slope populations. The optimum sample length 
depends on the size of the smallest topographic ele- 
ment of interest. The sample length was limited to 
approximately 0.75 km by the resolution of the pho- 
tographic plates and the aperture of the micropho- 
tometer. The effect of this sampling restriction on 
the slope-component measurements is best considered 
by examining schematic topographic profiles (fig. 9A) 
in which the individual slope elements are approxi- 
mately 0.75 km. A sample less than 0.75 km long 
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Figure 9. — Schematic topographic profiles: A, symmetrical 
topography with individual slope elements 0.75 km long; 
B, asymmetrical topography with slope elements 0.25 and 
0.50 km long. 

yields accurate average slope data, but one larger 
than 0.75 km would yield values lower than the actual 
slopes. A sample length of 1.5 km gives an average 
slope value of zero. In practice, alternate 0.75-km 
samples were recorded because of a limitation of the 
recording equipment (paper tape punch). This pro- 
cedure prevents accurate profile construction, but 
does not affect relative-roughness rankings. 

Topographic profile sampling presents problems 
similar to those inherent in function sampling. In 
general, the sample interval should be no longer than 
one-half the wavelength of the sampled function 
(Nyguist frequency) to avoid producing “aliased” 
data (Blackman and Tukey, 1958, p. 32). In figure 
9 B, the topographic profile is asymmetric ; the sample 
lengths here must be shorter than the shorter ele- 
ment (0.25 km) to insure adequate sampling. More 
complicated arrangements are obviously possible, 
both with terrestrial and lunar topography. 

In general, the number of measurements required 
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for adequate representation of a slope-component 
population depends on the roughness of the area. 
About 100 slope-component measurements generally 
suffice for apparently smooth relatively homogeneous 
areas; 200-300 are required for rough areas owing to 
the greater range of dispersion of slopes. 

Because of the limited availability and longer slope 
length (3-5 km) of the relative-relief data, their use- 
fulness is limited to the rougher terra terrain units 
(II-B, II-C, II-D). Although these data are useful 
in describing the relative roughness of the terra, the 
slope-component data are more widely applicable and 
will be discussed in more detail. 

SLOPE-FREQUENCY DISTRIBUTION FUNCTIONS 

Slope-frequency distribution functions are useful 
for comparison of slope-component populations be- 
cause the shape of the distribution curve depends on 
the frequency of slopes in each class interval. An 
early application of slope-frequency distributions to 
lunar terrain classification was made by McCauley 
(1964, 1965), who ranked the six terrain units of the 
1 :2, 000, 000 scale terrain map by comparing the char- 
acter of the typical calculated cumulative frequency 
distribution curve for each unit (fig. 10). 



Figure 10. — Cumulative slope-component frequency distribu- 
tion curves for the six terrain units (I through VI) shown 
on the 1 : 2,000,000-scale terrain map of the lunar equatorial 
zone. (After McCauley, 1964, fig. 9.) 

Statistics, such as the median (50 percentile value) 
and the 10-90 percentile range, are easily derived 
from these cumulative frequency curves and are use- 
ful for estimating relative-roughness relations. For 
example, it was shown during early systematic sam- 
pling of the equatorial zone that median slope values 
over 0.75-km slope lengths range from 5° for the 
heavily cratered terrae to less than 1° for typical 


lunar maria. The 10-90 percentile range, a measure 
of the dispersion of the slope values, is also a useful 
measure of roughness. In general, as the median 
slope values increase so also does the dispersion. For 
smooth maria, the dispersion is about 2°, and for the 
roughest lunar terrain sampled it is about 8°. Also, 
as the terrain becomes rougher, the area under the 
slope-frequency curve becomes greater as the curve 
shifts to the right with increasing median and dis- 
persion values, as shown in figure 10. 

Although slope-frequency data commonly are 
treated in the absolute sense (Strahler, 1954, p. 12; 
McCauley, 1964), certain advantages are realized by 
considering the slope-component direction, as well as 
the magnitude. In so doing the slope-component pop- 
ulations can be arranged in the form of an algebraic 
frequency distribution, rather than an absolute 
frequency distribution. 

The advantages gained by treating the slope- 
component data with respect to direction as well as 
magniture are threefold: (1) the algebraic frequency 
distributions generally approximate Gaussian, or 
normal, distributions and are, therefore, described by 
two statistics, the arithmetic mean and the standard 
deviation, (2) the frequency of topographic change 
can be estimated by measuring the number of times 
that the slope changes algebraic sign in a given pop- 
ulation (that is, the percentage of slope reversal 2 ), 
and (3) if the algebraic frequency distribution is 
non-Gaussian, the degree of asymmetry of the terrain 
can be estimated. The algebraic arithmetic mean and 
standard deviation, percentage of slope reversal, and 
the absolute arithmetic mean and median are given 
in table 2 for each of the areas sampled; the terrain 
units of each sample area and the location of the area 
are also listed. 

Algebraic slope-frequency distribution curves were 
derived for all the sample areas by using the <t(z) 
methods (Lindgren and McElrath, 1959, p. 253) for 
the standard normal distribution function. The fit is 
good for most of the sample areas. Notable excep- 
tions, however, are samples 20, 26, and 30 (fig. 8; 
table 2). These deviations appear to result from 
either inadequate sampling, censoring (as described 
on p. Gil), or superposition of significantly different 
terrain types. Percentile dispersion values for two 
typical sample areas (9, 39) are shown in figure 11 
to illustrate the Gaussian nature of algebraic alope- 
frequency distributions. For comparison, percentile 
dispersion values for non-Gaussian, absolute distri- 
butions are also shown. 

2 For example, if in a slope population sample of 100 measurements, the 
slope changes sign 10 times, the percentage of slope reversal is 20. 
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The value of algebraic slope-component distribu- 
tions for rapid visual appraisal of relative rough- 
ness of widely different lunar terrain types is dem- 
onstrated by examining typical mare and terra areas. 
The maria are characterized by a moderate density 
of craters in all size ranges and by widely spaced 
structures, such as ridges, domes, and rilles. Mare 
Cognitum (sample 9, the area photographed at 
greater resolution by Ranger VII) is an excellent 
example of this terrain type. The terrae include a 



SLOPE, IN DEGREES 

Figure 11. — Arithmetic probability plot of slope-com- 
ponent percentile values for typical samples, Mare 
Cognitum (sample 9) and Bode B (sample 39), 
showing the Gaussian character of the algebraic 
frequency distributions (solid line) and the non- 
Gaussian, absolute distributions (dashed line). (After 
Rowan and McCauley, 1966, fig. 9.) 

more complex variety of basic terrain types, from 
relatively smooth terra plains (generally resembling 
mare at 0.75 km resolution) to rougher units com- 
posed in great part of a continuum of fresh to ero- 
sionally degraded partly buried craters and irregular 
topographic high unrelated to craters. 


In figure 12, three different terrain types are 
contrasted by comparison of their algebraic slope- 
component frequency distributions. The relatively 
smooth Mare Cognitum sample (No. 9) is repre- 
sented by a leptokurtic curve at 0.75-km sample 
length and has 52 percent of the slope components 
confined within the ±1° range. This property is a 
common characteristic of many of the mare samples 
(for example, samples 8 and 12). In contrast, the 
terra (II-D) sample near Herschel (No. 51) exhibits 
a platykurtic curve. The frequency distribution here 
shows a greater range of slope-component values be- 
cause of the greater inherent roughness. In many 
areas, roughness is intermediate between that of the 
Mare Cognitum and Herschel samples. Sample 39, 
for example, exhibits a mesokurtic curve (Bode B, 
fig. 12). Thus, curves of this general type can be used 
to evaluate the relative roughness of widely different 
types of terrain, independent of geologic considera- 
tions. 

Algebraic slope-frequency distributions for dome, 
fields (samples 33, 47) are characteristically bimodal. 
The degree of separation between the modes depends 
primarily on the type of dome sampled. The flat sub- 
dued Hortensius domes show a small modal separa- 
tion, whereas the more complex Marius domes (not 
included in table 2) show a larger separation. 

Many of the mare samples contain several slope 
populations owing to the presence of features such 
as rilles and ridges. In the Gambart area (sample 

12) , hummocky terra remnants, small craters, rays, 
and low broad ridges are present on the mare (fig. 

13) , and the algebraic frequency curve is mesokurtic 
(fig. 14A ) . The degree to which the various popula- 
tions contribute to the roughness of the mare can be 
evaluated by excluding these populations individually 
from the total sample and observing the effects on 
the resulting slope-component frequency curves 
(14 B-F) . A moderately leptokurtic curve (fig. 14 B) 
results when the hummocky terra (fig. 13, unit 1) is 
deleted. Exclusion of the mare ridge (fig. 13, unit 2), 
the crater fields (fig. 13, unit 3), and the ray-covered 
areas (fig. 13, unit 4), as well as the hummocky terra 
(fig. 13, unit 1), shows that the character of the 
curves (fig. 14 C-E) is significantly affected by the 
superposed ray material, because it is at this point 
of exclusion (fig. 14 E) that the curve becomes mark- 
edly leptokurtic. Sampling of the mare area only (fig. 
13, unit 5; fig. 14 F), the only sample included in 
table 2, yields a highly leptokurtic curve. 

In contrast to algebraic slope-frequency distribu- 
tions, which are approximately Gaussian, the abso- 
lute slope-frequency distributions are positively 
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Table 2. Significant statistical slope-component parameters and the terrain units listed in decreasing order of percentage 
area covered for each of the 51 slope-component sample areas (slope lengths 0.75 km) 


IX, 


,, algebraic arithmetic mean; <r A] , algebraic standard deviation; X Ab , absolute arithmetic mean; Med Ah , absolute arithmetic mean] 


Sample area 

Sample 

Location 
(lat and long) 

Number of 
measurements 

*A> 

<T 

A1 

Percentage 

slope 

reversal 

X., 

Ab 

Med 

Ab 

Map units (table 1) 

Mare 

Fecund itatis. 

1 

1.2°-3.8°S., 

42.3°-44.7°E. 

900 

— 0.14° 

0.48° 

26.7 

0.40° 

0.32° 

I-A, IV-A, IV-B, I-B, III-C, 
IV-E, rays, III-A. 

Riphaeus 

Mountains. 

2 

8.1°-9.2°S., 

24.7°-26.7°W. 

582 

.38° 

.38° 

15.3 

.34° 

.29° 

I-A, IV-A, III-A-3, IV-E, 
III-A-1. 

Reiner __ _ 

_ 3 

5.3°-6.2°N., 

51.9°-54.9°W. 

1,428 

—.06° 

.62° 

15.9 

.47° 

.36° 

I-A, IV-A, III-A-3a, III-A-3b, 
and rays. 

Sinas A__ _ 

_ 4 

8.0°-9.5°N., 

32.1°-36.5°E. 

551 

—.20° 

.78° 

11.4 

.68° 

.64° 

I-A, IV-A, III-C, III-A-3, 
III-A-1, IV-E, and rays. 

Riphaeus 

Mountains. 

5 

8.3°-9.3°S., 

26.2°-27.5°W. 

402 

.56° 

.85° 

20.3 

.72° 

.43° 

• I-A, IV-A, IV-D, III-A-3, III-A-1, 
IV-E, and rays 

Lubbock S_. 

. 6 

1.3°-2.3°N. 

35.6°-39.4°E. 

509 

.00° 

.83° 

16.1 

.63° 

.51° 

I-A, III-C, II-B, II-C, 
and IV-G. 

Cayley 

7 

3.8°-4.9°N., 

14.7°-16.6°E. 

396 

—.51° 

.79° 

17.9 

.72° 

.56° 

II-A, II-B, III-A-3, III-A-1, 
III-A-2, III-C, III— B— 1 , III-B-2, 
III-B, IV-C-1, and rays. 

Sinus 

Medii. 

8 

0.75°-3.0°N., 

2.1°W-4.1°E. 

1,451 

.22° 

1.13° 

13.8 

.87° 

.69° 

I-A, IV-A, III-A-3, III-A-1, rays, 
II-C, IV-B, II-B, III-A, and 
IV-C-1. 

Mare 

Cognitum. 

9 

6.0°-10.0°S., 

8.9°-13.2°W. 

902 

.23° 

.99° 

17.7 

.79° 

.65° 

I-A, II-B, rays, III-A-1, III-A-3, 
III-C, IV-A, IV-D, IV-F, IV-E, 
III-B-3, III— B— 1 , and III-A-2. 

Mosting C__ 

10 

0.0°-2.0°S. 

8.15°-10.0°W. 

201 

—.21° 

1.10° 

16.4 

.89° 

.76° 

I-A, II-A, II-B, rays, and III-A-3. 

Sinus 

Aestuum. 

11 

8.0°-10.0°N.. 

10.0°-13.0°W. 

440 

.14° 

1.34° 

18.2 

1.08° 

.91° 

I-A, IV-A, rays, III-B-3, IV-C-2, 
IV-C-1, IV-H, and III-C 

Gambart 

12 

0.0°-3.0°N., 

10.0°-14.4°W. 

626 

.09° 

1.38° 

19.0 

1.11° 

.91° 

I-A, IV-A, III-B, rays, II-C, IV-E, 

III- A-3, III-A-1, III-A-2, 

IV- C-1, and IV-H, 

Theophilus GA 
SA 7-4. 

13 

5.3°-8.0°S. 

21.3°-28.4°E. 

360 

—.24° 

1.04° 

36.0 

.85° 

.71° 

III-A-3b and I-A. 

Theophilus GA 
SA 7-3. 

14 

5.3°-8.0°S., 

21.3°-28.4°E. 

678 

—.57° 

1.41° 

25.0 

1.12° 

.87° 

III-A-3b, II-B, and IV-H. 

Schroter G 
SA 6. 

15 

1.0°-3.1°N., 

7.0°-10.1°W. 

493 

.15° 

1.10° 

22.7 

.83° 

.62° 

I-A, IV-A, III-C, rays, IV-E, 

III- A-3, III-A-1, III-B-3, 

IV- B, and III-A. 

Schroter F 

16 

7.1°-10.0°N.,. 

3.0°-8.0°W. 

603 

.19° 

1.35° 

27.2 

.87° 

.57° 

I-A, IV-A, III-C, III-A, IV-B, 
II-B, III-B-3, III-B-1, II-A, 
and rays. 

Dionysius 

17 

0.3°-1.3°N., 

16.5°-17.3°E. 

361 

—.04° 

.121° 

33.2 

.83° 

.54° 

II-B and II-A. 

Flammarion 

18 

2.3°-4.0°S., 

2.7°-4.2°W. 

226 

—.27° 

1.55° 

22.0 

1.09° 

.89° 

III-B-2a, III-C, III-B, IV-C-1, 
III-A, and rays. 

Kunowsky 

19 

0.0°-2.5°N., 

31.6°-34.6°W. 

969 

.07° 

1.81° 

14.5 

1.15° 

.95° 

I-A, rays, II-C, III-C, III-A-3, 
III-A-I, and III-A-2. 

Reaumur __ 

20 

1.7°-3.3°S., 

0.1°W-1.1°E. 

162 

.03° 

1.80° 

19.0 

1.06° 

.48° 

III-B-2, III-B, II-A, and IV-C-1. 

Lansberg D 

21 

2.5°-5.6°S., 

30.0°-32.5°W. 

185 

.15° 

1.60° 

17.8 

1.27° 

.97° 

I-A, IV-A, III-A-3, III-A-1, 
II-B, IV-B, rays, III-A-2, and 
IV-E. 

Arago E 

22 

9.25°-8.0°N., 

21.3°-25.5°E. 

822 

—.17° 

1.56° 

23.1 

1.21° 

.98° 

I-A, IV-A, III-C, III-A-1, III-A-3, 
IV-B, and rays. 

Schroter F 

23 

7.1°-10.0°N., 

3.0°-8.0°W. 

188 

—.19° 

1.56° 

25.0 

1.17° 

.84° 

I-A, IV-C-1, II-B, III-A, III-C, 
III-B, and IV-D. 
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Table 2. — Significant statistical slope-component parameters and the terrain units listed in decreasing order of percentage 
area covered for each of the 51 slope-component sample areas (slope lengths 0.75 km) — Continued 


Sample area 

Sample 

Location 
(lat and long) 

Number of 
measurements 

X A1 

<7 

A1 

Percentage 

slope 

reversal 

X Ab 

Med 

Ab 

Map units (table 1) 

Hortensius (2) 
SA 1. 

24 

7.1°-8.1°N., 

26.6°-29.4°W. 

349 

.17° 

2.01° 

1.52 

15.2° 

1.23° 

I-A, rays, and IV-A. 

Euclides F 

_ 25 

3.0°-6.2°S., 

31.6°-35.0°W. 

1,026 

.53° 

1.80° 

22.6 

1.44° 

1.08° 

I-A,IV-A, rays, III-A-3, III-A-1, 
III-C, and IV-E. 

Hipparchus 

, 26 

4.7°-6.2°S., 

3.0°-6.9°E. 

404 

.30° 

1.38° 

24.2 

1.03° 

.75° 

\ 

III-B-2, III-B, IV-C-1, IV-B, rays, 
III-B-1, III-B-3, III-A-3, 

III-A-1, III-A-2, and III-A. 

Ptolemaeus 

. 27 

7.6°-10.0°S., 

0.5°E-3.6°W. 

805 

—.06° 

1.69° 

27.0 

1.14° 

.83° 

III-B-2a, III-B, III-B-2, III-B-1, 
III-B-3, III-A-3, III-A-1, 
III-A-2, and IV-H. 

Hyginus Rille 
SA 1. 

28 

7.0°-10.0°N., 

4.9°-9.9°E. 

1,042 

.10° 

1.93° 

21.8 

1.41° 

1.10° 

II-A, III-C, III-A, IV-C-1, 
III-A-3, II-C, and III-A-1. 

Hyginus Rille 
SA 2. 

29 

7.0°-10.0°N., 

4.9°-9.9°E. 

248 

—.25° 

1.93° 

27.4 

1.49° 

1.22° 

II-B, IV-C-1, III-A-3, III-A-1, 
and III-A. 

Fra Mauro___ _ 

. 30 

5.2°-6.9°S., 

15.5°-18.0°W 

282 

.09° 

1.84° 

21.2 

1.51° 

1.35° 

II-B, II-D, III-B-2a, IV-C-1, 
III-A-3, III-A-1, rays, III-B-1, 
III-C, and III-B-3. 

Lansberg DB___ 

. 31 

3.1°-4.1°S., 

35.0°-31.6°W. 

569 

- — .23° 

1.98° 

25.1 

1.49° 

1.16° 

I-A, IV-A, III-A-3, III-A-1, rays, 
and III-C. 

Theophilus GA 
SA 6. 

32 

5.3°-8.0°S Ii 

21.3°-28.4°E. 

560 

.00° 

1.64° 

27.1 

1.14° 

.78° 

I-A, IV-H, IV-F, and III-A. 

Fauth G 

. 33 

4.0°-6.0°N., 

14.0°-16.0°W. 

148 

.12° 

1.79° 

27.7 

1.41° 

1.17° 

II-C, I-B, IV-B, IV-C-1, rays, 
I-A, III-C, and III-A-3b. 

Turner 

. 34 

0.0°-1.0°S., 

10.0°-14.0°W. 

272 

—.21° 

1.46° 

21.0 

1.18° 

.98° 

I-A, rays, IV-A, II-B, II-C, 
IV-D, IV-E, III-A-1, III-A-3, 
III-C, and III-A. 

Pallas D-_ 

. 35 

1.5°-3.2°N., 
2.6°— 4.5°W. 

224 

.30° 

2.14° 

22.7 

1.62° 

1.40° 

I-A, II-A, IV-A, IV-D, II-C, 

II- C, IV-E, III-C, III-A-3, 

III- A-1, III-B-3, III-B-1, III-A, 
and III-A-2. 

Silberschlag P_ - 

. 36 

6.0°-7.5°N., 

10.0°-12.0°E. 

145 

—.11° 

2.50° 

26.2 

1.94° 

1.46° 

II-B, III-C, III-A, III-A-3, 
III-A-1, and rays. 

Glyden 

. 37 

5.0°-6.0°S., 

0.0°-1.2°E. 

79 

.24° 

2.53° 

29.0 

1.89° 

1.43° 

II-C, IV-A-3, II-D, III-A-2, 
III-A-1, and rays. 

Schroter F 

- 38 

7.1°-10.0°N., 

3.0°-8.0°W. 

590 

— .57° 

2.31° 

30.8 

1.83° 

1.51° 

II-B, III-A-3, III-A-1, III-A-2. 
IV-C-1, III-B, IV-H, III-B-3, 
III-B-1. IV-G, III-C, and III-A. 

Bode B 

. 39 

7.0°-9.0°N., 

1.0°—4.5°W. 

449 

.06° 

2.81° 

23.8 

2.20° 

1.82° 

II-B, I-A, III-A-3, III-A-1, 
III-A-2, II-A, III-C, III-B-3, 

III- B-1, rays, III-A, III-B, 

IV- C-2, and IV-C-1. 

Schroter F . . . 

. 40 

7.1°-10.0°N., 

3.0°-8.0°W. 

430 

—.23° 

2.42° 

34.9 

1.84° 

1.46° 

II-B, III-B, III-A-3, III-A-1, 
III-C, IV-G, III-B-3, III-B-1, 
and IV-C-2. 

Schroter F 

. 41 

7.1°-10.0°N., 

3.(T-8.0 o W. 

218 

—.44° 

2.63° 

32.1 

2.03° 

1.66° 

II-B and II-C. 

Theophilus GA 
SA 7-2. 

42 

5.3°-8.0°S., 

21.3°-28.4°E. 

77 

.06° 

2.53° 

35.0 

1.98° 

1.53° 

II-B, III-A-3b, and III-C. 

Schroter A9 

. 43 

5.1 0 -6.6°N., 

7.9°-9.5°W. 

594 

—.46° 

2.80° 

32.2 

1.92° 

1.37° 

II-B, I-A, IV-A, III-C, III-B-3, 
and III-B-1. 

Fauth 

. 44 

4.0°-6.0°N., 

16.0°-20.0°W 

4,096 

20.° 

2.96° 

24.4 

2.20° 

1.76° 

i rv oo 

II-C, I-B, III-A-3b, IV-G, III-A, 
IV-H, and IV-C-1. 

TT T» TT rr TTT A TTT n TTT A O 


Hyginus Rille 45 7.0°— 10.0°N., 864 .35° 2.87° 23.5 2.29° 1.98 II B, II C, III A, III C, III A 3, 

SA 3. 4.9°-9.9°E. III-A-1, IV-H, and IV-G. 
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Table 2. Significant statistical slope-component parameters and the terrain units listed in decreasing ovdev of percentage 
area covered for each of the 51 slope-component sample areas (slope lengths 0,75 km ) — Continued , 


Sample area 

Sample 

Location 
(lat and long) 

Number of 
measurements 

*A1 

o 

A1 

Percentage 

slope 

reversal 

X.. 

Ab 

Med 

Ab 

Map units (table 1) 

Theophilus GA 
SA 1. 

46 

5.3°-8.0°S., 

21.3°-28.4°E 

225 

—.35° 

2.83° 

32.9 

2.12° 

1.63° 

III-A-3b and IV-C-1. 

Hortensius (2) 
SA 2. 

47 

7.1°-8.1°N., 

26.6°-29.4°W. 

77 

—.04° 

3.68° 

29.9 

2.92° 

2.53° 

IV-B, II-C, IV-H, IV-A, and rays. 

Schroter G 

SA 3. 

48 . 

1.0°-3.1°N., 

3.0°-8.0°W. 

401 

—.58° 

3.74° 

29.4 

2.20° 

1.35° 

II-B, III-C, IV-C-1, III-B-3, 
III— B— 1, and rays. 

Schroter F 

49 

7.1°-10.0°N., 

3.0°-8.0°W. 

270 

.51° 

3.63° 

38.1° 

2.87° 

2.32° 

II-C, II-B, I-A, III-A-3, and 
III-A-1. 

Theopilus GA 

50 

5.3°-8.0°S., 

21.3°-28.4°E. 

382 

—.12° 

4.29° 

22.0 

3.54° 

3.23° 

II-D, III-B-3, IV-H, and III-B-1. 

Herschel 

51 

5.0°-6.6°S., 

4.0°-2.7°W., 

and 

1.3°W-1.2°E. 

290 

0 

00 

o 

T 

4.06° 

29.4 

3.21° 

2.69° 

III-A-3,. II-D, II-C, III-B, 
III-A-3, III-A-1, III-A-2, and 
III-C. 
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Figure 13. — Telescopic photograph (Pease No. 192, Mount 
Wilson Observatory) of the Gambart area (sample 12), 
with the nature and extent of different slope-component 
populations outlined. 1, hummocky terra; 2, mare ridge; 3, 
crater field; 4, ray-covered area; 5, mare. (After Rowan 
and McCauley, 1966, fig. 11.) 

skewed and are, therefore, non-Gaussian. Curves fit- 
ted to the absolute frequency distributions are “re- 
versed J-shaped curves,” characteristic of folded nor- 
mal frequency distributions (Leone and others, 
1961). Accurate estimates of central tendency and 
dispersion values from non-Gaussian frequency dis- 
tributions require transformation to the standard 
normal distribution. To avoid this laborious task, 
which could be facilitated by computer processing, 
algebraic rather than absolute slope-frequency dis- 
tribution curves were used as the principal terrain 
rating values in this study. 

SLOPE-COMPONENT STATISTICS 
In describing frequency distribution functions, 
two general types of statistics are commonly used: 
central tendency _and dispersion values. Both the 
arithmetic mean, X, and the median are central ten- 
dency measures, whereas the standard deviation, cr, 
is the most widely used expression to indicate dis- 
persion of population distributions. Somewhat less 
commonly used values are the mode, the class interval 
with the highest frequency, and the 10-90 percentile 
dispersion. Of occasional value are the skewness, 
which is an expression of the departure from the 
normal frequency distribution, and the kurtosis, 
which measures peakedness of the distribution. 

To assess the value of each of these statistics as 
measures of relative roughness, all were calculated 
and evaluated for each of the 51 slope-component 


sample areas. Listed in table 2 are those statistics 
that appear most useful; the subscripts, Ab and ai, 
indicate values taken from absolute and arithmetic 
slope-component frequency distributions, respec- 
tively. 

Evaluation of all the standard statistical values 
of the slope-component samples suggests that those 
most useful as relative-roughness parameters are the 
absolute arithmetic mean, Xu>, the algebraic standard 
deviation, ct-ai, and the percentage of slope reversal; 
of somewhat more limited use is the algebraic arith- 
metic mean, Xu. The meaning of these statis- 
tics in terms of relative roughness is demonstrated 
by the schematic topographic profiles in figure 15, 
which show progressively decreasing wavelengths, 
W, and increasing amplitudes, A, and, consequently, 
progressively rougher terrain. To completely de- 
scribe these profiles, only representative values of W 
and A are necessary. If the sampling procedure could 
be controlled so that the sample length was always 
W/2 and in phase, W and A would effectively de- 
scribe the surface. Controlled sampling of this type 
was, of course, not possible in this study, and more- 
over, terrain profiles are rarely simple single func- 
tions. Thus, the slope-component samples are only 
statistical representatives of the real roughness of 
the terrain units previously described on the 
1:1, 000, 000-scale map. 

In the schematic profiles, the topographic texture 
or frequency could be estimated by the percentage of 
changes of direction of the slope from positive to 
negative (that is, the percentage of slope reversal). 
The percentage of slope reversal would, at a constant 
sampling cell length and rate, be greater in the high- 
frequency (rough) terrain (fig. 15C) than in the low- 
frequency (smooth) areas (fig. 15A, B) . Slope re- 
versal values are particularly useful because they 
provide a means of distinguishing between slope- 
component populations which have similar arithme- 
tic means and standard deviations but different 
topographic textures. 

Slope and relative-relief measurements obtained at 
the same sample length are, of course, simply related 
and can be of equal value in expressing relative 
roughness, but slope-component statistics were the 
principal ranking tool because of the ease by which 
they are generated from photoclinometry. As pre- 
viously mentioned, the arithmetic mean, X, is a mea- 
sure of the central tendency and is, therefore, the 
average value of the slope-component frequency dis- 
tribution. The absolute arithmetic mean, Xu>, of the 
schematic profiles (fig. 15) increases as the slopes in- 
crease, but the algebraic arithmetic mean, Xu, will 
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SLOPE, IN DEGREES 


Figure 14. — Algebraic slope-component frequency distribution curves for the Gambart sample area (No. 12) with various 
slope-component populations shown in figure 13 deleted: (A) total sample area; (B) area 1 deleted; (C) areas 1 and 2 
deleted; (D) areas 1 and 3 deleted; (E) areas 1 and 4 deleted; and (F) areas 1, 2, 3, and 4 deleted. (After Rowan and 
McCauley, 1966, fig. 12.) 

be zero unless the profiles are tilted or the topogra- 
phy is markedly asymmetric. The algebraic standard 
deviation, o- A i, a measure of dispersion of the alge- 
braic slope-frequency distribution, increases pro- 
gressively as the slope increases in the profiles. In 
terrain more complex and diversified than that rep- 
resented by these schematic profiles, no single sta- 
tistical value obtained at a constant sample cell 
length appears to adequately express relative-rough- 
ness differences, and combinations of various param- 
eters are required. Although the absolute arithmetic 
mean, Xu,, describes the average slope, the algebraic 
standard deviation, cr A i, is needed to express the di- 
versity of the slopes, and the texture of the terrain is 
best described by the percentage of slope reversal. 
The usefulness of slope-component statistics as 



c 


Figure 15. — Schematic topographic profiles showing the 
increasing roughness as the wavelength (W) decreases 
and the amplitude (A) increases from A to C. 
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relative-roughness indicators is demonstrated by 
analysis of the Schroter F area, which contains both 
terra and mare. Six terrain samples were delineated 
(fig. 16) by inspection prior to photoclinometric sam- 
pling. Comparison of the Xui, cr A i, and percentage of 



Figure 16. — The Schroter F area — samples 16, 23, 38, 40, 41, 
and 49 — showing the boundaries of the slope-component 
samples and the morphologic types plotted on a telescopic 
photograph (Pease No. 124, Mount Wilson Observatory). 
(After Rowan and McCauley, 1966, fig. 7.) 

slope reversal shows that sample areas 16 and 23 are 
very similar, as are areas 40, 41, and 49, and area 
49 is obviously the roughest. If prior qualitative 
analysis of the terrain had not been made, these dis- 
tinctions would have been masked (tables 3, total 
area), because more than one slope-component popu- 
lation would have been included in the sample. 

To gain a better appreciation of the significance 
of central tendency slope-component values over 
varying slope lengths, slope studies of three moun- 
tainous terrestrial areas were made. Selection was 
based in part on the availability of good topographic 
maps at increasing scales and on general random- 
ness of slope direction, as expressed by a predomi- 
nantly dendritic drainage pattern. The Sierra 
Nevada— Tuolumne Meadows area in Yosemite Park, 
Calif., the San Juan Mountain area north of Du- 
rango, Colo., and the Mogollon Rim region between 
McNary and Payson, Ariz., were studied at three 
scales: 1:1,000,000 (Aeronautical Chart and Infor- 


mation Center map) , 1 :250,000 (Army Map Service 
map), and 1:62,500 (U.S. Geological Survey map). 
In each region, slope components were measured at 
slope lengths of 1.0, 0.25, and 0.06 km. Only alternate 
slope-component segments were used to simulate as 
closely as possible the sampling procedure used in 
the microphotometer. The data were then aranged 
in cumulative frequency diagrams (fig. 17) for com- 
parison with lunar slope-component data used in 
the early phase of this study. At the 1:1,000,000 
scale (slope length = 1 km), the median slope- 
component values of these samples of rough moun- 
tainous terrestrial terrain range from 2.5° to 4.0° 
(fig. 17A). At larger scales, however, the median 
increases. For example, at a slope length of 0.25 
km the median for the for the Sierra Nevada sample 
is 6.0° (fig. 17B), and at 0.06 km it is 8.5° (fig. 17C) . 
Absolute median slope-component values for the 
lunar terrain samples are listed in table 2 for com- 
parison with these terrestrial values. 

Three relationships that bear on the general 
problem of terrain roughness emerge from these 
terrestrial studies. First, the measurable roughness 
clearly is dependent on the length or the interval 
over which slope or any other topographic property 
is measured. The limiting slope length used for sam- 
pling purposes is clearly a function of the resolution 
of the photography employed or the resolution of 
the topographic map from which slope data are 
derived. Second, as resolution increases, the rough- 
ness, as expressed by central tendency values or 
other slope statistics, also increases. Third, and most 
significant, low statistical slope parameters over 
long slope lengths do not mean that these surfaces 
are necessarily smooth at larger scales. The surfaces 
may contain many locally smooth areas at larger 
scale, but the statistical roughness for the entire 
area sampled will increase rather than decrease, at 
a rate characteristic of the terrain unit itself. Each 
terrain unit may also be categorized by its own 
distinctive rate of increase of roughness with in- 
creasing resolution. A further discussion of this 
problem is contained under “Extrapolation tech- 
niques.” 
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Figure 17. — Cumulative slope-frequency distribution 
curves for rough terrestrial terrain measured at 
varying scales: A, 1:1,000,000, 1.0-km slope length; 
B, 1:250,000, 0.25-km slope length; and C, 1:62,500, 
0.06-km slope length. (After McCauley, 1964, fig. 
11 .) 


RELATIVE-ROUGHNESS PARAMETERS 

To evaluate the relative roughness of the major 
1 : 1,000,000-scale lunar terrain units at earth-based 
resolutions, the absolute arithmetic mean, Xib, the 
algebraic standard deviation, cr A i, and the percent- 
age of slope reversal calculated for the 51 sample 
areas were plotted agarnst each other. The most use- 
ful plot appears to be Xu> versus o- A1 (fig. 18) . Mare 
and terra divisions clearly occupy discrete zones or 
areas on the plot, although there is some overlap; 
crater components span the range of roughness from 
mare to rugged terra (units II-C and II-D). Of 
particular importance is the clearly different rough- 
ness of unconstructed smooth mare relative to mare 
complicated by ridges,' rilles, and craters. The ranges 
of characteristic slope-component and relative-rough- 
ness parameters, X Ab , <r A i, and percentage of slope 
reversal, are given in table 4. 

According to Leone, Nelson, and Nottingham 
(1961), Xu, for “half -normal distributions” is re- 
lated to o r A i by 

Xu, = yj-^- • o- A i. (1) 

Although figure 18 clearly shows that the two pa- 
rameters are closely related, equation 1 does not hold 
for the data plotted in figure 18 ; the Xu, values are 
too low for given o- A i values. An empirical expression 
relating these parameters for photoclinometrically 
derived lunar slope components was derived from 
figure 18: 

X Ab = 0.79 or^ ; (2) 

Xu, can be estimated from equation 2 and need not 
be derived independently, and <x A i can be used by it- 
self as an index of relative roughness and as an ex- 
pression for the rate of change of roughness as a 
function of scale when plotted against the sample 
length. 

_ Plots of the percentage of slope reversal versus 
X Ab and cr A i exhibit a large amount of scatter ; thus, 
slope reversal may be only indirectly related to these 
statistical values. Slope reversal may, therefore, be 
a fundamental terrain classification parameter in the 
same sense that amplitude and frequency are critical 
in the definition of a periodic function. 

As previously mentioned, Xu proved to be of little 
general use. Where topography is symmetrical in the 
east-west direction ami there is no censoring because 
of shadowing effects, Xai will be close to zero. A sig- 
nificant departure from zero in the absence of the 
above effects, however, is a measure of the regional 
slope, which is almost impossible to detect by quali- 
tative study of the lunar plates. Thus, X A i has local 
significance. In areas where regional slopes are 
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Table 4 . — Limiting relative-roughness values for major terrain units in the lunar equatorial zone 




Slopo-component statistics 



Relative relief (meters) 
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Percentage of 
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Figure 18. — Plot of the algebraic standard deviation, <r A i, versus the absolute arithmetic mean, Y Ab , for the 51 slope-compo- 

hent sample areas showing the range of the major lunar terrain types. 
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greater than 1 °^_ Xu has been adjusted to zero in 
table 2, so that Xu> is not affected. 

Because of the previously described difficulties in 
obtaining dependable slope-component measurements 
in relatively rough areas, relative-relief statistics 
were used as supplemental means of describing rela- 
tive roughness in the lunar terra. Approximately 
1,200 measurements, for sample lengths ranging 
from 3 to 5 km, were statistically analyzed for the 
three roughest terra units, II-B, II-C, and II-D, and 
rougher mare, I-A. The measurements were fur- 
nished by the Air Force Chart and Information Cen- 
ter, St. Louis, Mo. Table 4 gives the most useful 
relative-relief statistical parameters for the terrain 
types where data were available, and the frequency 
distributions are shown in figure 19. These distribu- 
tions are all bimodal, suggesting that the measure- 
ments were made on at least two distinct terrain 
populations. The terra areas in the equatorial zone 
are commonly a composite of several types of terrain 
consisting of (1) rugged hills and intervening valleys 
of structural origin, (2) hummocky to braided ter- 
rain related presumably to the Imbrium event, and 
(3) superimposed old and young craters. Relative- 
relief data clearly discriminate hummocky terra, 
II-B, from the more rugged sculptured terra, II-C 
and II-D. Effective detailed quantitative classifica- 
tion is, however, difficult owing to this morphologic 
heterogeneity. During regional mapping of these 
terrain types, qualitative estimates of relative rough- 
ness had to be relied upon. 

Different empirical combinations of these basic 
relative-roughness parameters can be used to classify 
lunar terrain in terms of suitability for mobility and 
landing. For example, in considering a landing site, 
roughness classification might be weighted in favor 
of the frequency data because numerous abrupt slope 
reversals in any area would present a severe hazard. 
On the other hand, if mobility were the prime con- 
sideration, the classification might be weighted in 
favor of the standard deviation, a measure of ter- 
rain diversity ; a high standard deviation would indi- 
cate the presence of numerous areas of extreme slope 
that could not be traversed. 

EXTRAPOLATION TECHNIQUES 

The rate of change of roughness, as defined by one 
or more terrain parameters plotted as a function of 
the size of the sample cell (A L), appears to be of 
value for limited extrapolation of terrain character- 
istics below the resolution of particular photographs. 
It may be distinctive for various units, and for to- 
pography characterized by increasingly smaller cra- 
ters of the same forms, may remain constant over a 





RELATIVE RELIEF, IN METERS 

c 

Figure 19. — Relative-relief frequency histograms for A, 
hummocky terra, II-B; B, hummocky terra with moder- 
ate relief, II-C; and C, sculptured terra with moderate 
to high relief, II-D. N is the number of measurements. 
Measurements made over 3-5 km horizontal distances by 
the U.S. Air Force Chart and Information Center, St. 
Louis, Mo., using the shadow technique. (After Rowan 
and McCauley, 1966, fig. 14.) 
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wide range of resolutions. Other areas photographed 
under progressively higher resolution will show in- 
flection points in the curve as new types of terrain 
features, such as large blocks, are resolved. 

This extrapolation technique was used to predict 
the mean slope of typical maria at the 1-meter scale 
from earth-based observations prior to the flight of 
Ranger VII (McCauley, 1964) . Figure 20 shows that 


ers (1965) as 1-meter A L. They are very similar to 
the parameters and curves for terra at the 1-km 
scale or resolution (table 2; figs. 12, 18). Statistics 
derived from the original photometric scans and 
those from the map are not, however, strictly com- 
parable. 

Although few other data points are currently 
available, it appears that predictive techniques might 



Figure 20. — Extrapolation of median slope to the 1-meter scale from telescopic data for mare terrain. The 
segment WX represents extrapolation prior to receipt of Ranger VII data. Curves for each distinctive 
terrain unit should show a different slope. (After McCauley, 1964, fig. 14; and Rowan and McCauley, 
1966, fig. 15.) 


the median slope-component value at the 1-meter 
scale (about 3.6°), derived from the original photo- 
metric scans of the last Ranger VII frame (P 3 -979) , 
falls near a straight-line extrapolation of median 
slopes over 1.0 km and longer slope lengths (Rowan 
and McCauley, 1966). 

Figure 21 gives the statistical relative-roughness 
parameters Xu> and <x A i and slope frequency curve 
derived from analysis of profiles of the map prepared 
of the last Ranger VII frame by Shoemaker and oth- 


be further developed by examining the functional 
relations between A L and other statistical values 
such as crater density, relative-relief data, slope re- 
versal, and power spectral data for a variety of lunar 
terrain types. Wood and Snell (1959) investigated 
the use of a ranking method of correlation for six 
terrain factors from 200 samples in the United 
States. Average large-scale slope could be predicted 
from small-scale relief maps with an error of 4 per- 
cent for the smoother parts of the United States. 
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Figure 21. — Algebraic slope-component frequency distribution curve and statistical parameters for Ranger VII photograph 
P3-979 of an area in Mare Cognitum. (After Rowan and McCauley, 1966, fig. 16.) 
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Using curvature rather than slope measurements, 
personnel at the Bendix Systems Division (1965) 
also demonstrated the relations between A L and 
roughness. Once the effect of resolution (AL) on the 
statistical parameters of particular lunar terrain 
units is established, large-scale surface properties 
may be predictable with a fair degree of accuracy 
if the terrain unit can be identified in lower resolu- 
tion photography. Lunar Orbiter monoscopic 1-3- 
meter resolution photographs of many terrain types 
are now available, and 60-90-meter resolution pho- 
tographs of nearly all the lunar surface are also 
available. Although a large sampling program will 
be required, nomograms can be constructed for most 
terrain units to permit the prediction of relative 
roughness in areas not covered by 1-3-meter reso- 
lution photographs. The data that can be obtained 
by terrain extrapolation techniques could be very 
useful for planning future planetary exploration 
missions; the techniques, therefore, deserve addi- 
tional study. 

APPLICATIONS TO SPACECRAFT TARGETING 

As previously mentioned, the primary objectives 
of the lunar terrain analysis program were to clas- 
sify and rank terrain in the equatorial zone accord- 
ing to relative roughness so that the results could 
be applied to site-selection studies for the Ranger, 
Surveyor, Lunar Orbiter, and Apollo missions. Dur- 
ing the latter stages of the program a rationale was 
developed for the selection of photographic targets 
and landing sites. Essential to this rationale were 
data on the relative roughness of the major terrain 
units in the equatorial zone and an appreciation for 
the predictable greater roughness with increasing 
resolution. When statistical relative-roughness data 
for a potential landing site can be extrapolated from 
low-resolution data, the relative probability of a 
successful landing can be predicted. In addition, data 
from landed spacecraft can be extrapolated to a po- 
tential landing site, if similar terrain types occur 
there. It is impossible, however, to determine in ad- 
vance whether unresolved hazards exist at a particu- 
lar' point. 

RANGER 

The successful Ranger experiments, VII, VIII, and 
IX, obtained nested, increasing-resolution television 
images of two widely separated maria and the floor 
of a modified terra crater. Results of the terrain 
analysis studies were used in the selection of the 
Ranger VII and Ranger VIII targets. The selected 
aim point for Ranger VII was long 20° W. and lat 
11° S., in the unnamed mare between Mare Nubium 
and Oceanus Procellarum (this area was later 
named Mare Cognitum in recognition of the success 


of Ranger VII) . This first view of the moon at the 
1-meter scale provided by Ranger VII indicated that 
its surface, although heavily cratered, was probably 
smooth enough for spacecraft landings. Ranger 
VIII was successfully deployed to photograph a small 
area in Mare Tranquillitatis (long 24° E., lat 3° N.) . 
Ranger IX was sent to the floor of the terra crater 
Alphonsus mainly to gather data on an area in which 
transient activity had been reported and also to 
photograph terra basin materials. 

The Ranger spacecraft, in addition to providing 
views of very small areas on the moon at a resolution 
of about a meter, were important because they pro- 
vided continuous variable-scale imagery that bridged 
the gap between earth-based telescopic resolution 
and the type of large-scale data needed for landing- 
site evaluations. The large-scale high-resolution data 
returned by the three Ranger spacecraft suggested 
that areas of the moon were sufficiently smooth at 
the 1-meter scale for landings. The feasibility of both 
the Surveyor and Apollo programs was thus partially 
demonstrated. The principal task of the Lunar Or- 
biter and Surveyor missions was to extend the high- 
resolution Ranger coverage and thereby provide the 
basis for certifying an array of equatorial sites 
suitable for manned landings. 

SURVEYOR 

Prior to the launch of Surveyor I, 65 candidate 
landing sites were selected (McCauley and others, 
1966) on the basis of telescopic observations and 
photographs, terrain and geologic maps, the statis- 
tical relative-roughness data reported herein, and 
the results of analyses of the Ranger VII and VIII 
photographs. 

The rationale for the selection of these landing 
sites was based on various spacecraft engineering 
constraints, scientific merits, and terrain roughness. 
The 65 sites were also rated by a partially quantita- 
tive technique to facilitate final value judgments 
whereby the engineering, scientific, and site suita- 
bility considerations were traded off against one 
another. 

LUNAR ORBITER 

The Lunar Orbiter spacecraft were designed to 
photograph potential landing areas for Surveyor and 
Apollo. Each of the five spacecraft was equipped 
with a dual-lens camera capable of obtaining mono- 
scopic 1-meter and stereoscopic 8-meter resolution 
photographs from a perilune altitude of approxi- 
mately 50 km. 

Planning the photographic tasks of the Lunar Or- 
biter missions was a complicated undertaking be- 
cause of lighting, exposure, and readout problems. 
The specific requirements of the Surveyor and Apollo 
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programs for the best possible photography of po- 
tential landing areas near the equator further re- 
stricted targeting capability. At the time of initial 
planning for Orbiter I, the only photography avail- 
able with resolution greater than about 0.5 km was 
that obtained by Rangers VII, VIII, and IX and Sur- 
veyor I. Although of great scientific and engineering 
value, the high-resolution photographic coverage of 
these experiments was very limited in areal extent. 
Moreover, the area photographed by Rangers VII 
and IX were, for the most part, outside of the equa- 
torial zone. Therefore, during mission planning ac- 
tivities extensive use was made of the 1:1, 000, 000- 
scale terrain map (Rowan, 1966; Holm and others, 
1965) and the lunar geologic maps then available. 

Most of the major terrain types in the equatorial 
zone were photographed during the first Orbiter mis- 
sion, including the dark maria, all four types of 
terra, and modified and well-formed craters. Detailed 
terrain maps of each of nine prime sites to be pho- 
tographed by Lunar Orbiter I were prepared just 
prior to the mission and used during the mission for 
locating more precisely the desired photographic 
coverage. (See U.S. Geol. Survey, 1966.) Each suc- 
ceeding mission plan incorporated the results of 
evaluation of the preceding missions, but as new 
sites were selected, the 1 : 1,000,000-scale terrain and 
geologic maps continued to be used in the site- 
selection process. 

Because the high-resolution camera of Orbiter I 
failed early in the mission, the two successive Orbit- 
ers were also devoted almost exclusively to the land- 
ing-site search. Elimination of large tracts of the 
equatorial belt based on terrain consideration facili- 
tated effective completion of the site-search role 
prior to Orbiter IV. The last two missions could be 
devoted, therefore, to gathering new scientific data, 
first by means of Orbiter IV's planetwide coverage 
of most of the lunar surface at an average resolution 
of 80 meters and later by Orbiter V’s detailed 
2-meter photography of sites of particular scientific 
interest. 

SUMMARY 

Results of the lunar terrain analysis program 
were as follows: (1) Lunar equatorial terrain was 
classified and terrain maps were prepared at 
1:2,000,000 and 1:1,000,000 scales, (2) relative- 
roughness indicators were derived through statis- 
tical analysis of relative-relief data and some 150,000 
slope-component measurements at 0.75-km slope 
lengths made during this study, (3) the usefulness 
of algebraic slope-frequency distributions was dem- 


onstrated at telescopic resolutions, and (4) relative 
roughness of mare surfaces at small slope lengths 
was extrapolated from roughness data for larger 
slope length, and thereby the utility of the tech- 
nique for future landing and trafficability studies 
was demonstrated. In addition, the photoclinometric 
technique was refined and automated for slope meas- 
urements from telescopic photographs so that enough 
measurements were available for valid statistical 
analysis. Perhaps the most promising aspect of these 
results is the application of the extrapolation tech- 
niques to quantitative description of terrain and 
geologic units. 

These results have been applied to site-selection 
studies for the photographic Ranger and Lunar 
Orbiter missions and for Surveyor and Apollo land- 
ings. Much of the rationale developed in the site- 
selection studies could be applied to the exploration 
of other terrestrial planetary surfaces. 
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4 Figure 22 - — Representative lunar terrain units in southern 
part of Marius Hills area near crater Reiner {lat 8° N., 
ion g 53° W. ; Lick Observatory plate ECD 75), 

Unit I~A. Smooth mare plains without resolvable relief 
<X) ; low crater density. 

Unit IV-A, Asymmetrical mare ridge (2); smoother 
symmetrical mare ridge (3). 

Unit IV-B. Low symmetrical dome with side slopes of 
several degrees and relief of about 100 meters (4) ; 
steep-sided dome with side slopes up to 10° and 
relief of about 200 meters (5). 

Unit IV-D. Broad flat- topped mare ridge or plateau 
bounded by gentle escarpments with slopes of several 
degrees (6), 
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FIGURE 24,— Representative lunar terrain units in in ter mare 
area northwest of Fra Mauro (lat 1,5 c S, p long 20 Q W.; 
Mount Wilson Observatory plate Pease 124 C) : 

Unit II-B. Hummocky to gently rolling terrain typical 
of much of region surrounding Mare Imbrium (1) ; 
slopes about 5° on flanks of hummocks or linear hills. 
Unit III-A. Sharp unmodilied crater (2), 
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< 4 Figure 23.— Representative lunar terrain units in southern 
Sinus Medii area (lat 0 D N,, long 2° £,; Lick Observa- 
tory plate ECD 63), 

Unit I- A, Smooth mare plains (1), 

Unit I-R. Dark mare plains; few craters; no visible 
relief (2). 

Unit 11-A. Smooth terra plains; higher density in 
small craters than in mare units; planar to gently 
rolling surface (3), 

Unit II-C. Blocky to lineated terra; smooth at' fine 
scale but slopes on individual blocks about 10°-20 fl 
(*)■ 

Unit III-A. Fresh unmodified crater; sharp rim crest; 
surrounded by rays; interior slopes up to 30°; ex- 
terior slopes on rim about 15* (5), 

Unit IV-C-I. Linear rilles, bounded by steep relatively 
straight scarps with slopes about 15° -20° (6) ; fiat 
floors. 
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4 Figure 25* — Representative lunar terrain units in area 
between craters Ptolemaeus and Hipparchus (lat 7.5° S., 
long 3* E. ; Lick Observatory plate ECD 65) : 

Unit II- A. Terra plains (1), 

Unit IV— 1-1. Chain craters (2) ; coalescing craters 
with low smooth rims located along a linear structure. 


0 50 100 KILOMETERS 

1 J 


M Figure 26*— Representative lunar terrain units in area 
northeast of crater Ptolemaeus (lat 7 D S. f long 6 b W*; 
Lick Observatory plate ECD 65) : 

Unit II- A* Terra plains (1) ; gently rolling in places. 
Unit II-C- Hummocky to blocky, lineated terra (2)* 


Figure 27* — Representative lunar terrain units in Apen- ► 
nine Mountains, south of Mare Imbrium (lat 25° N. p long 
3° E.; Lick Observatory plate ECD 35), 

Unit I-B* Dark mare (1)* 

Unit II-D. Rough high blocky terra (2)* Side slopes 
about 20“ and relief of 3,000-4,000 meters. 

Unit IY-C-2. Sinuous rille (3); steep interior slopes, 
flat floor. 
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A Figure 28, — Representative lunar terrain units in AristiU 
lus area in Mare Imbrium (lat 35° N n long 0 s E; Lick 
Observatory plate ECD 35) : 

Unit III-A-la. Steeply sloping, terraced inner crater 
wall (1) ; slopes up to 30°. 

Unit III-A-2a, Smooth crater floor (2), 

Unit I II — A— 2c, Crater central peaks (3) ; rough blocky 
slopes up to 20 s , 

Unit III-A-3a* Rough terraced crater rim (4) ; slopes 
up to 15* ; extends outward for about one- half of the 
crater diameter* 

Unit III-A-3b, Radial ridge terrain on crater rim; 
extends outward for about one crater diameter (5) ; 
markedly smoother than unit III-A-3a. 
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A Figure 29. — Representative lunar terrain units in Alba- 
tegnius area (lat 12° S., long 4* E,; Lick Observatory 
plate ECD 18) : 

Unit IIL-B-1* Subdued or modified crater wall (1); 
slopes about 10* ; smooth terraces. 

Unit III-B-2c, Subdued crater central peaks (2) ; 
smoother and of lower relief than peaks in unmodi- 
fied large craters. 

Unit IXI-B-S. Rolling, gently terraced to hummocky 
terrain on crater rim (3); traceable outward for one- 
quarter to one-balf of the rim diameter. 
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